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ABSTRACT 

We present high-speed, multi-colour optical photometry of the anomalous X-ray pulsar 
4U 0142+61, obtained with ULTRACAM on the 4.2-m William Herschel Telescope. 
We detect 4U0142+61 at magnitudes of i' = 23.7 + 0.1, q' = 27 .2 + 0.2 and v! > 25.8, 
consistent with the magnitudes found bv lHulleman et alTl)2004|) and hence confirming 
their discovery of both a spectral break in t he optical and a l ack of long-term optical 
variability. We also confirm the discovery of iKern fc Martini pOO^ that 4U 0142+61 
shows optical pulsations with an identical period (~ 8.7s) to the X-ray pulsations. The 
rms pulsed fraction in our data is 29 + 8%, 5-7 times greater than the 0.2-8 keV X-ray 
rms pulsed fraction. The optical and X-ray pulse profiles show similar morphologies 
and appear to be approximately in phase with each other, the former lagging the 
latter by only 0.04 ± 0.02 cycles. In conjunction with the constraints imposed by X- 
ray observations, the results presented here favour a magnetar interpretation for the 
anomalous X-ray pulsars. 
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1 INTRODUCTION 

More than 100 X-ray pulsars are currently known. The vast 
majority of these are found in low-mass and high-mass X- 
ray binaries (LMXBs and HMXBs), and are hence pow- 
ered by accretion onto a rotating, magnetised neutron star. 
There exists a small group of 8 X-ray pulsars, however, 
that exhibit properties very much at variance with those 
of the accreting pulsars in X-ray binaries. These so-called 
Anomalous X-ray Pulsars (AXPs) all have ~ 5 — 12 s spin 
periods which decrease steadily with time, soft (and rela- 
tively low-luminosity) X-ray spectra, no radio emission, and 
tend to be associated with supernova remnants in the galac- 
tic plane. Most importantly, the AXPs show no evidence 
of a binary companion. Fo r a recent review of AXPs, see 
IWoods fc Thompson! <2004h . 

This latter fact prompted a variety of models based on 
isolated neutron stars and white dwarfs (see the review by 
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llsrael et alJl2002D . but these run into difficulty on energetic 
grounds: the loss of rotational energy, which powers radio 
pulsars like the Crab, is orders of magnitude too small to 
power the observed X-ray luminosity of the AXPs. An addi- 
tional energy source is therefore required, for which two com- 
peting models seem to have emerged: accretion from a fossil 
disc or ultra-strong magnetic fields. In the former scenario, 
an isolated neutron star accretes from a fossil disc, such 
as might be produced through fall-back of material after a 
supernova explosion or left over from a common-envelope 
phase which destroyed the companion star. In the latter 
scenario, AXPs are "magnetars", isolated neutron stars with 
enormous (B ~ 10 14 — 10 15 G) magnetic fields. It is the decay 
of the magnetic field which heats the neutron star surface, 
causing it to emit thermal radiation in the X-rays. Non- 
thermal emission is then produced by particles accelerated 
in the magnetosphere by the Alfven waves from small-scale 
fractu res on the neutron star surface llThompson fc Duncanl 
1996) or Compt onization of thermal ph otons by magneto- 
spheric currents llThompson et alll2002l) . 
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The magnetar model has begun to dominate the litera- 
ture in recent years. There are sound theoretical reasons for 
why this is so, as the magnetar model now appears to be able 
to explain the observational properties of several categories 
of supposedly young neutron stars that are not powered 
by rotation, including the AXPs and the Soft Gamma Re- 
peaters (SGRs). Such unification is sup ported by the recent 
discovery of SGR -like bursts in AXPs jGavriil et al.ll200l 
iKaspi et al1l2003l) . which suggests there might be an evolu- 
tiona ry link between AXPs and SGRs (see lMereehetti et alJ 
120021 and references therein) . 

But what other observational evidence is there to sup- 
port the magnetar model of AXPs? Belief in the magnetar 
model rests partly on the failure of the accretion model to 
explain the faintness of the optical/infrared counterparts, 
which sets strong limi ts on the size of an accretion disc (e.g. 
iHulleman et alJ l2000l) . and the fact that the pulsed frac- 
tion of optical light is significantly greater than it is in X- 
rays, ruling out reproc essing of X-rays in a disc as its origin 
<Kern fc Martini 12001 but see lErtan fc ChendEooi) . Both 
of these optical constraints have been obtained via obser- 
vations of the brightest known AXP, 4U 0142+61. In this 
paper we report on new high-speed, multi-colour optical ob- 
servations of this object, obtained with the a im of confirming 
the hi gh optical pulsed fraction observed bv lKern fc Martini 
(2002). 



2 OBSERVATIONS AND DATA REDUCTION 

The observations of 4U 0142+61 presented in this paper 
were obtained wit h ULTRACAM dDhillon fc MarsblboOll 
iBeard et alJ 120021) at the Cassegrain focus of the 4.2-m 
William Herschel Telescope (WHT) on La Palma. ULTRA- 
CAM is a CCD camera designed to provide imaging photom- 
etry at high temporal resolution in three different colours 
simultaneously. The instrument provides a 5 arcminute field 
on its three 1024 x 1024 E2V 47-20 CCDs (i.e. 0.3 arcsec- 
onds/pixel). Incident light is first collimated and then split 
into three different beams using a pair of dichroic beamsplit- 
ters. For the observations presented here, one beam was ded- 
icated to the SDSS u (3543A) filter, another to the SDSS 
g (4770A) filter and the third to the SDSS (7625A) i' fil- 
ter. Because ULTRACAM employs frame-transfer chips, the 
dead-time between exposures is negligible: we used ULTRA- 
CAM in its two-windowed mode, each of 100 x 200 pixels, 
resulting in an exposure time of 0.48 s and a dead-time of 
0.025 s. A total of 30 618 and 31304 frames of 4U 0142+61 
were obtained on the nights of 2002 September 10 and 12, 
respectively, with each frame time-stamped to a relative ac- 
curacy of better than 50 [is using a dedicated GPS system. 
Both sets of data were obtained in photometric conditions, 
with no moon and i'-band seeing of 0.75 and 0.65 arcseconds 
on 10/09/02 and 12/09/02, respectively. 

A portion of the summed i' image from the night of 

1 The absolute timing accuracy of ULTRACAM was verified 
with contemporaneous observations of the Crab pulsar. Our ob- 
served time of optic al pulse maxim um was found to agree with 
the ephemeris of lLvne et alJ 120051) to better than 1 millisecond 
(the quoted error in the Crab pulsar ephemeris during September 
2002). 



Table 1. Updated ephemeris for 4U 0142+61 spanning the opti- 
cal observations descri bed in sect ion |21 ba sed on the monitoring 
campaign described in iGavriil fc Kaspil l2002Tl . BMJD refers to 
the Barycentric-corrected Modified Julian Date on the Barycen- 
tric Dynamic al Timescale (TDB). TO A refers to the pulse time 
of arrival (see lGavriil fc Kaspil <2002l) for details). The errors on 
the last two digits of each parameter are given in parentheses. 



BMJD range 51610.636 - 53 401.184 

TOA arrival points 79 

v (Hz) 0.11509507445(18) 

v (10~ 14 Hz s" 1 ) -2.66478(36) 

v (10- 24 Hz s- 2 ) 5.08(23) 

Epoch (BMJD) 52 506.9748874274228 

rms residual (cycles) 0.031 



12/09/02 is shown in figure The vertical streaks are due 
to light from bright stars falling on the active area of the chip 
above the CCD windows. The data we obtained on 10/09/02 
(not shown in figure suffer from streaks passing through 
the position of 4U 0142+61, increasing the background noise 
level significantly. As a result, we rotated the Cassegrain 
rotator in advance of our observations on 12/09/02 so that 
no streaks passed through 4U 0142+61. For this reason, the 
data obtained on 12/09/02 are of a much higher quality 
than the data obtained on 10/09/02. Note that the vertical 
streaking problem has since been rectified in ULTRACAM 
by provision of an adjustable focal-plane mask which blocks 
the light from bright star s (and the sky) above the CCD 
windows (IStevensor]l2004D . 

The data were reduced using the ULTRACAM pipeline 
software. All frames were first debiased and then flat-fielded, 
the latter using the median of twilight sky frames taken with 
the telescope spiralling. We then extracted light curves of 
4U 0142+61 using two different techniques: 

2.1 Technique (i) 

iKern fc Martini (|2002) obtained their light curve of 
4U 0142+61 by synchronising the CCD clocks in their cam- 
era to the X-ray spin period of 4U 0142+61, resulting in the 
accumulation of 10 on-chip phase bins. This has the advan- 
tage of reducing detector noise, but the potential disadvan- 
tage that a period must be assumed before the data have 
been taken and, if the period is wrong, the true pulse profile 
is unrecoverable. 

To mimic the IKern fc Martini (I2002T) technique, we 
assumed a spin period for 4U 0142+61 on 12/09/02 of 
8.688473130 s, which was calculated from the updated X- 
ray ephemeris given in table 1. Note that this ephemeris 
spans our WHT observations and is hence more reliable for 
our p urposes than using the ephemeris of FGavriil fc Kaspil 
( 2002) adopted by IKern fc Martini ll2002l) . Each ULTRA- 
CAM data frame was then added to one of 10 evenly-spaced 
phase bins covering the spin cycle of 4U 0142+61, resulting 
in 10 high signa l-to-noise dat a frames. An optimal photom- 
etry algorithm (|Navlorl ll998') was then used to extract the 
counts from 4U 0142+61 and a bright comparison star 24 
arcseconds to the east of the AXP (see figure 0, the latter 
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Figure 1. Left: Summed i' image from the night of 12/09/02, with a total exposure time of 15046 s (~ 4 h). Star A is 4U 0142+61 and 
star B is the comparison/reference star (sec Hullc man et alJ |2Q04) for coordinates and magnitudes). The orientation arrows represent 
10 arcseconds on the sky. For clarity, only a portion of the two ULTRACAM windows is shown. Note that there is no gap between the 
windows, but a faint discontinuity between them can be seen running down the centre of the image, due to the fact that each window is 
read out via a separate channel. Right: Higher contrast plots of the field around 4U 0142+61 (star A), showing the summed i' (top) and 
g 1 images (bottom) from the night of 12/09/02. The box in the left-hand image shows the portion of the field shown (at the same scale) 
in the right-hand images. 



acting as the reference for the profile fits and transparency- 
variation correction. The position of 4U 0142+61 relative to 
the comparison star was determined from a sum of all the 
images, and this offset was then held fixed during the re- 
duction so as to avoid aperture centroiding problems. The 
sky level was determined from a clipped mean of the counts 
in an annulus surrounding the target stars, and subtracted 
from the object counts. 



2.2 Technique (ii) 

The second approach we took to light curve extraction was 
identical to that described above, except that we omitted the 
phase-binning step and simply performed optimal photom- 
etry on the 61922 individual ULTRACAM data frames. In 
other words, we made no assumption about the spin period 
of 4U 0142+61. 



3 RESULTS 
3.1 Magnitudes 

We were unable to detect 4U 0142+61 in it', at a detection 
limit of u > 25.8. We did, however, clearly detect it in g 
and i' on both nights at magnitude s of q' = 27.2 ± 0.2 and 
i' = 23.7 ± 0.1, as shown in figure^ H ulleman et alJ (120041) 



measured g' ~ 26.9 and i' ~ 23.7 (where we have con- 
verted their BVRI Johnson- Morgan- Cousins magnitudes 
to SPSS m agnitudes using the transformation equations of 
ISmith et ai]l2002ft . indicating that 4U 0142+61 was approx- 
imately the same magnitude during our observations. 

3.2 Pulse profiles 

The two data reduction techniques described in section [5] 
result in two different pulse profiles for 4U 0142+61. 

3.2.1 Technique (i) 

The first technique produced the pulse profiles shown in 
the top panel of figure [5] As expected, the light curve of 
12/09/02 is of a significantly higher quality than that of 
10/09/02, but both show approximately the same mor phol- 
ogy as the optical pulse profile presented bv lKern fc Martini 
(2002), exhibiting a broad (arguably double-humped) struc- 
ture with peaks around phases 0.65 and 1.15 and a mini- 
mum around phase 0.35. These phases are different to the 
corres ponding phases in the pulse profile of iKern fc Martini 
(2002) , but this is t o be expected given that, as discussed 
bv IKern fc Martini i2002T) . their optical observations were 
obtained outside the span of the ephemer is they used and 
the s ource does exhibit some timing noise iGavriil fc Kaspil 
2002). Our timing solution, on the other hand, is based on 
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figure[5]can be measured using a peak-to-trough pulsed frac- 
tion, h p t, denned as follows: 

, Fmax Fmin /~\ 

V = f — Tf~' (1) 

* max I ± min 

where Fmax and F m i n are the maximum and minimum 
flux in the pulse profile, respectively. We find a value of 
hpt — 58±f6% on f 2/09/0 2, higher than the puls ed fraction 
of h pt = 27±8% derived bv lKern fc Martini §002), although 
the difference between the two values is only marginally sig- 
nificant (31 ± 18%, i.e. < 2a). There are a number of factors 
which might contribute to a higher optical pulsed fraction 
in our data: 
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Figure 2. Top: Pulse profiles of 4U 0142+61 in the i'-band on 
12/09/02 (solid line), obtained using technique (i) (section 12.11 . 
The dotted line shows the poorer quality light curve we obtained 
on the night of 10/09/02, demonstrating the repeatability of the 
pulse profile. Each pulse profile was first corrected for trans- 
parency variations using the comparison star (star B in figure^, 
although the correction made only a negligible difference to the 
light curves. The pulse profiles were then normalised by divid- 
ing by the mean number of counts. Note that the formal error 
bars on these pulse profiles were unreliable (most probably due 
to the vertical streaks shown in figure^, and hence the error bars 
shown have been calculated from the scatter in the light curve ex- 
tracted using technique (ii) ("section ^, 21 . Bottom: Averaged X-ray 
pulse profile of 4U 0142+61 in the 2-10 keV en ergy band, which 
is an u pdated version of the profile presented in lQavriil fc Kasnil 
12002). Note that it is not possible to estimate the X-ray pulsed 
fraction from this profile as the background level (i.e. the mini- 
mum flux) in the X-ray pulse profile is unrelated to the pulsar - 
see section|l]for details). For this reason, no scale is given on the 
ordinate. 



the updated ephemeris for 4U 0142+61 presented in table 1, 
which spans our WHT observations and is hence reliable. 

There is some similarity in the morphologies of the opti- 
cal pulse profile shown in the top panel of figure[5]and the 2- 
10 keV X-ray pulse profile shown below it, where the latter is 
an up dated version of the data presented bv lGavriil fc KaspH 
(2002). Both profiles share a similar broad/double- humped 
morphology. Moreover, since the X-ray light curve shown in 
figure [5] has also been phased using the ephemeris given in 
table 1, it can be seen that the optical and X-ray pulse pro- 
files are approximately in phase with each other. To quan- 
tify this, the optical pulse profile was cross-correlated with 
the X-ray pulse profile. The resulting peak in the cross- 
correlation function was fitted with a parabola to derive a 
shift of 0.04 + 0.02 cycles (i.e. 0.35 + 0.17 s), where a positive 
phase shift implies that the optical pulse profile lags the X- 
ray pulse profile. This result is only marginally significant 
(at the 2a level), due to the low signal-to-noise and time 
resolution of the optical data, and additional data will be 
required in order to confirm that the phase shift is signifi- 
cantly different from zero (discounting the unlikely situation 
in which the time delay is approximately equal to some mul- 
tiple of the spin period) . 

The modulation amplitude of the pulses presented in 



• The pulsed fraction mea sured from our data re fers to 
the i' band, whereas that of iKern fc Martini l|200 2*l is for 
white light (4000-10000A). If the optical pulsed fraction 
varies with wavelength, this could be the source of the dis- 
crepancy. Note that our g' data were too faint to extract a 
pulse profile from, unfortunately, so we are not in a position 
to test this explanation. 

• Even a small error in the assumed period on which the 
data is phase-binned can result in a smearing of the pulse 
profile and hence a reduction in the measured pulsed frac- 
tion. We have simulated this effect and fin d that to reduce 
our pu lsed fraction to the level observed bv lKern fc Martini 
(2002), the period must be in error by greater than ~ 
0.003 s. This is three orders of magnitude greater than the 
timing accuracy achie ved by the instrumentation used by 
IKern fc Martini ll2002t) and hence an error in the period used 
to phase bin the data is an unlikely source of the discrepant 
pulsed fractions. 

• The higher pulsed fraction in our data might be due ei- 
ther to a decrease in the unpulsed optical c omponent or an 
increa se in the pulsed optical component. iHulleman et alJ 
(2004) found no evidence for long-term i?-band variability 
in 4U 0142+61, down to a 2-a limit of 0.09 magnitudes, 
and our magnitude estimates (section 13.11 app ear to sup- 
port t his conclusion. Note, however, that IHulleman et alJ 
(2004) did find long-term variability of ~ 0.5 magnitude 
in their isT-band observations of 4U 0142+61, which they 
tentatively attributed to the occurence of SGR-like bursts 
in 4U 0142+61. X-ray observations of this source have not 
shown such bursts, but this might be due to their (expected) 
low amplitude and the efficiency of the X-ray monitoring. 

• The peak-to-trough pulsed fraction defined in equa- 
tion 1 effectively adds any noise present in the light curve 
to the true pulsed fraction, thereby tending to increase the 
resulting measurement. A more robust estimate is given by 
the root-mean-square (rms) pulsed fraction, h rms , defined 
as follows: 



hrms — _ 



1=1 



(2) 



where n is the number of phase bins per cycle, yi is the 
number of counts in the i th phase bin, a t is the error on 
yi and y is the mean number of counts in the cycle. As 
expected, measuring the optical pulsed fraction in this way 
gives a lower value of h rms = 29 ± 8% . This is much closer 
to the value derived bv lKern fc Martini <2002l) . but it should 
be stressed that these authors measured the peak-to-trough 



ULTRACAM observations of the AXP 4U 01^2+61 5 




0.1 1 

frequency (Hz) 

Figure 3. Lomb-Scargle periodograms of 4U 0142+61 in the i'- 
band on 10/09/02 (top panel) and 12/09/02 (bottom panel), ob- 
tained using the light curves from technique (ii) (section 12.21 , 
The dotted line shows the predicted X-ray pulse frequency of 
0.11509502130 Hz on 12/09/02, calculated from the ephemeris 
given in table 1. 



pulsed fraction (equation 1), not the rms pulsed fraction 
(equation 2). 

• The higher pulsed fraction in our data could be due 
to some systematic problem with the sky subtraction. We 
consider this to be unlikely, however, as one would not 
th en expect our magnit ude estimates to agree with those 
of iHulleman et alJ (120041) . 



3.2.2 Technique (ii) 

The second data reduction technique (section 12. 21 can be 
used to provide a check on the reliability of the optical pulse 
profile shown in figure |2] To do this, it is first necessary to 
fold the extracted light curve on the pulse period. Rather 
than do this by adopting the X-ray ephemeris given in ta- 
ble 1, as we did in figure |2] we can instead determine the 
pulse period directly from our optical data using a peri- 
odogram and then fold the data on this period. 

Figure shows the Lomb-Scargle periodograms 
JPress fc Rvbickilll989l) for the 30 618 and 31 304 points in 
the i' light curves obtained on 10/09/02 and 12/09/02, re- 
spectively. The g' light curves were unfortunately too noisy 
to perform such an analysis. The light curves were first 
corrected for transparency variations and then detrended 
by subtracting their mean level. The highest peak in the 
resulting periodogram of 12/09/02 occurs at a period of 
8.687 ± 0.002 s, where the error is given by the width (a) of 
a Gaussian fit to the peak in the periodogram. This period 
is consistent with the X-ray pulse period given in table 1. 
Although noisier, an equivalent peak is also present in the 
periodogram of 10/09/02, with a period of 8.688 ± 0.002 s, 
thereby confirming that we have indeed detected the X-ray 
pulsation of 4U 0142+61 in the optical. We further tested 
the robustness of our period detection by constructing 10000 
randomised light curves from the original light curves by 
randomly re-ordering the y-axis points. Only 0.12% of the re- 
sulting 10000 periodograms for the 12/09/02 dataset showed 
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Figure 4. Pulse profiles of 4U 0142+61 in the i'-band on 
12/09/02, obtained using technique (ii) (section 12.21 and folding 
the resulting data on the optically-determined period of 8.687 s 
(solid line) and on the X-ray period of 8.688473130 s (dotted line). 
The data were first corrected for transparency variations using the 
comparison star (star B in figure 0. The pulse profiles were then 
normalised by dividing by the mean number of counts. Note that 
the formal error bars on these pulse profiles were unreliable (most 
probably due to the vertical streaks shown in figure^, and hence 
the error bars shown have been calculated from the scatter in the 
unfolded light curve. 



a higher peak at 8.687 s, and only 0.38% showed a higher 
peak at 8.688 s in the 10/09/02 dataset. 

Folding the i' light curve of 12/09/02 on the derived 
optical pulse period of 8.687 s gives the pulse profile shown 
in figure 2] The same data folded on the X-ray pulse period 
of 8.688473130 s is shown for comparison. Note that the 
phasing of both profiles can be directly compared to that in 
figure |1J as all of the data were folded using the zero point 
given in table 1. 

As one would expect, the i'-data folded on the X-ray 
pulse period (dotted line in figure is in excellent agree- 
ment with that presented in figure^ in terms of morphology, 
phasing and pulsed fraction (h pt — 50 ± 20%). The i'-data 
folded on the optically-determined pulse period (solid line 
in figure shares approximately the same phase of pulse 
maximum and pulsed fraction (h p t = 56 ± 16%), but the 
morphology is slightly different. In particular, the shape and 
phase of pulse minimum is very different, This is to be ex- 
pected, however, given that the data have been folded on 
the optically-derived period of 8.687 s, which is much less 
accurate than the X-ray period due to the lower quality of 
the optical data. 



4 DISCUSSION AND CONCLUSIONS 

Using two different data reduction techniques we have shown 
that the optical light from 4U 0142+61 pulsa tes on the X-ray 
period , thereby confirming the discovery of iKern fc Martini 
(2002). The morphologies of the 2-10 keV and i' pulse 
profiles are quite similar, both exhibiting a broad/double- 
humped structure. The optical lags the X-rays by only 
0.04 ± 0.02 cycles (0.35 ± 0.17s), i.e. there is no strong evi- 
dence for a phase shift between the two pulse profiles. The 
most reliable value we have derived for the optical pulsed 
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fraction is h rms — 29 ± 8%, as this is the more robust 
rms figure (as opposed to the peak-to-trough value) and has 
been obtained by folding the best dataset, that of 12/09/02, 
on the accurately known X-ray period (see figure . The 
X-ray pulsed fraction of 4U 0142+61 cannot be measured 
from the X-ray pulse profile shown in figure |5] unfortu- 
nately, as the background level (i.e. the minimum flux) in 
these data is unrelated to the pulsar. This is because the 
X-ray data were obtained with the Proportional Counter 
Array (PC A) on the Rossi X-ray Timing Explorer (RXTE), 
which has approximately a 1° field of view a nd no imag- 
ing ca pability. Instead, we turn to the work of iPatel et alJ 
(2003), who reported X-ray pulsed fractions for 4U 0142+61 
of h rms = 4.6+0.5% between 0.2-1.3 keV, h rms = 4.1+0.4% 
betw een 1.3-3.0 keV an d h rms — 5.6 + 1.0% between 3.0-8.0 
keV; IPatel et al] i2003T) also quote corresponding peak-to- 
trough pulsed fractions of h p t = 8.4+1.6%, h p t = 7.4 + 1.2% 
and hpt = 11.7 ± 3.2%. These values demonstrate that the 
optical rms pulsed fraction we have measured is 5-7 times 
greater than the X-ray rms pulsed fr actions, consistent with 
the factor of 5-10 times derived by iKern fc Martini (|2002T l 
from their peak-to-trough pulsed fraction measurements. 

We have measured g' and i' m agnitudes consistent wit h 
the BVRI magnitudes found by iHulleman et al.l l)2004l) . 
supporting their finding that, although variable in the in- 
frared and X-rays, 4U 0142+61 does not appear to show 
long-term variability in the optical part of the spectrum. 
Moreover, the fact that w e have used the B and V mag- 
nitudes of IHulleman et all ]2004l) to derive a g magnitude 
consistent with our own c onfirms that the spect ral break 
between B and V found bv IHulleman et all (120041) is real. 

The optical observations presented in this paper 
the refore lend addi ti onal w e ight to the arguments given 
bv IHulleman et al] feOOfft . IKern fc Martini (|2002ft and 
IHulleman et al] J2004I) that the AXP's are best explained 
by the magnetar model, mainly thanks to the failure of most 
alternative models to explain the observatio ns (a notable ex- 
cepti on is the disc-star dynamo gap model o flErtan fc Chenel 
2004) . In particu lar, fall-back accretion disc models (e.g. 
Perna et al] 1200011 fail because the optical flux is assumed 
to be due to reprocessing of the X-ray flux in the disc and 
therefore would not be expected to show either an optical 
pulsed fract ion significantly in ex cess of the X-ray pulsed 
fraction Csee lKern fc Martinl20 02) or a non-thermal spectral 
energy distribution in the optical. In addition, such repro- 
cessing might also be expected to result in the optical pulses 
lagging the X-ray pulses in phase by an amount depending 
on the light-travel time to the reprocessing structure, the re- 
processing timescale within it and its location with respect 
to the X-ray source and Earth. We have shown that there 
is no strong evidence for a phase shift between the optical 
and X-ray pulses, lending further evidence in support of the 
magnetar model. 



William Herschel Telescope is operated on the island of La 
Palma by the Isaac Newton Group in the Spanish Obser- 
vatorio del Roque de los Muchachos of the Instituto de As- 
trofisica de Canarias. 
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